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Abstract 

We investigate the algebraic structure of a recently proposed integrable t — J 
model with impurities. Three forms of the Bethe ansatz equations are presented 
corresponding to the three choices for the grading. We prove that the Bethe ansatz 
states are highest weight vectors of the underlying gl(2\l) supersymmetry algebra. 
By acting with the gl(2\l) generators we construct a complete set of states for the 
model. 
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1. Introduction 



The recent interest in the study of impurities coupled to both Fermi and Luttinger liquids 
has lead to the introduction of new theoretical approaches in order to study such systems. 
One approach is to adopt the prescription of the Quantum Inverse Scattering Method 
(QISM) which permits impurities to be incorporated into a model in such a way that 
integrability is not violated. The origins of this method trace back to the work of Andrei 
and Johannesson [1] who constucted an integrable version of a spin 1/2 Heisenberg chain 
containing an arbitrary spin S impurity (see also 0, [|). 

The long studied t — J model describes strongly correlated electrons interacting on a 
lattice comprising nearest neighbour hopping and spin exchange interactions. The model 
is defined on a restricted Hilbert space such that double occupancy of electrons at a 
single site is forbidden. On a one- dimensional lattice with periodic boundary conditions 
there exists a special choice of t, J such that the model can be solved exactly J|, [|, |6|]. 
Subsequently, it was shown that at this point the model acquires supersymmetry with 
respect to the Lie superalgebra g/(2|l) 0. The supersymmetry reflects the fact that 
the model possesses both bosonic and fermionic degrees of freedom. The existence of 
solvability and supersymmetry can be understood by the fact that the model can be 
derived in the context of the QISM || |9] which also proves integrability. Application of the 
algebraic Bethe ansatz leads to three different sets of Bethe ansatz equations corresponding 
to different choices of the grading. Furthermore, theoretical techniques based on the 
underlying supersymmetry expose significant consequences for the multiplet structure of 
the model. Specifically, it was shown in || that the Bethe states are highest weight 
states of the gl(2\l) superalgebra and that a complete set of states could be obtained by 
exploiting this supersymmetry. 



Adopting the QISM, Bediirftig et. al. [10] introduced an integrable periodic t — J 



model with impurity. In this model, the nature of the impurity is such that doubly 
occupied states can occur at those sites at which the impurities are situated. The origin 
of the impurity comes from changing the representation of the gl(2\l) generators at some 
lattice sites from the fundamental three-dimensional representation to the one parameter 
family of typical four dimensional representations which were introduced in |Tl[ to derive 
the supersymmetric U model. The resulting impurity model is significant in that the 
parameters from the four dimensional representations can be varied continuously without 
losing integrability. The impurities couple to both spin and charge degrees of freedom 
and these parameters describe the strength of the coupling of the impurities to the host 
chain. Limiting cases describe effective supersymmetric t — J models with additional sites 
or decoupled impurities in a t — J chain with twisted boundary conditions. Using the 
method of the algebraic Bethe ansatz, a solution for a particular choice of grading was 
derived. Also, a comprehensive analysis of thermodynamic properties such as ground state 
structure, elementary excitations, low-temperature specific heat and transport properties 
was given in |TU] . 

In this paper, we undertake an investigation into the detailed algebraic aspects of this 
impurity model. First, we present all three forms of the algebraic Bethe ansatz solutions 
corresponding to the three different choices of grading; viz. FBB, BBF, BFB. The case of 
FBB is that given in [IC|. A significant point here is that only recently have new Bethe 
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ansatz methods been proposed in order to solve the cases of the other gradings because of a 
lack of a suitable (unique) reference state. Rather one is forced to work with a subspace of 
reference states. This approach was developed in the works of Abad and Rios for an SU(3) 

TM to find a Bethe ansatz 



impurity spin chain [12|, [L3| and has already been adopted in 
solution of the supersymmetric U model starting from a ferromagnetic space of states and 
in the works |T^, [T7| regarding another integrable t — J model with impurities. At 
this point we mention that other impurity t — J models have been studied. In the work of 
Bares impurities were introduced into the model by way of inhomogeneities in the transfer 
matrix of the system [18|. Kondo (magnetic) impurities for arbitrary spin S have been 
studied by Schlottmann and Zvyagin [fLlj . These works treat cases of bulk impurities 
with periodic boundary conditions. It is also possible to consider open chains where the 
impurities are situated on the boundary. In these cases the impurity interactions assume 
very different forms than the case of bulk impurities and examples of these models have 
been considered in works of Wang et al [20| and Zhou et al |H[ . 

We also show that the Bethe ansatz states are highest weight states with respect to the 
underlying g>Z(2|l) supersymmetry algebra. By acting with the gl(2\l) lowering operators 
on the Bethe states we obtain additional eigenvectors. Then by a combinatorial argument 
we claim that in fact we recover a complete set of states for the model. Aside from the 
(non- impurity) integrable t — J model, completeness of the Bethe ansatz solution has been 
studied for example for the cases of the Heisenberg model |23|, the Hubbard model 
24 1 and the g/(2|2) supersymmetric extended Hubbard model [p5| |. Completeness for the 



integrable spin S impurity Heisenberg model was discussed by Kirillov [ffifl . Our results 
give the first example of completeness for a supersymmetric, integrable impurity model. 

The paper is organized as follows. In section 2 we recall the construction of the model 
through the QISM and also fix notation. In section 3 we discuss the algebraic Bethe ansatz 
solutions of this model. The highest weight property of the Bethe states is demonstrated 
in section 4 and the completeness of the Bethe states is shown in section 5. 



2. Quantum Inverse Scattering Method 



Recall that the Lie superalgebra g/(2|l) has generators J=1 satisfying the commuta- 
tion relations 

[£*., E{] = 6*Ei - (-i)(ra+0)(M+M) ( $j.E* (1) 
where the Z2-grading on the indices is determined by 



for i = 1,2, 



[i] = 1 for i = 3. 
This induces a Z2-grading on the gl(2\l) generators through 



+ (mod2). 



The vector module V has basis {V}f =1 with action defined by 



E)v 



5 k 3 v\ 



(2) 
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Associated with this space there is a solution R(u) G End(V (g) V) of the Yang-Baxter 
equation 

R l2 {u - v)R l3 {u)R 23 {v) = R 23 {v)R l3 {u)R 12 {u - v) (3) 
on the space V ®V ®V which is given by 

R{u) = -lj2e)®ei(-l)W. (4) 

U i,j 

We remark that eq. (|3|) is acting on a supersymmetric space so the multiplication of tensor 
products is governed by the relation 

(a®6)(c®d) = (-l) [b][c] ac®bd (5) 

for homogeneous operators b, c. 

The solution given by eq. (||) allows us to construct a universal L-operator which reads 

L(u) =1®!-^— $>}®£f(-L) b1 (6) 



i, j 



for any complex parameter A. This operator gives us a solution of the Yang-Baxter 
equation of the form 

R 12 (u - v)L 13 (u)L 23 (v) = L 23 (v)L l3 (u)R 12 {u - v) 

on the space V (g) V ® gl(2\l) which follows from the commutation relations eq. (|IJ). We 
use the one-parameter family of four-dimensional representations of gl(2\l) |TTJ acting on 
the module W with basis {v a , v b , v c , v d }. The explicit matrix representatives are given by 



E\ = 


\b) (c\ , 


El = 


\c) (b\ , 


E\ = 


- |c> (c\ - \d) (d\ , 


El = 


- \b) (b\ - \d) (d\ , 


El = 


\fa \a) (b \ + \/ a + 1 c) (d\ , 


E\ = 


\fa \b) (a\ + \/a+ 1 \d) (c , 


E\ = 


— \fa \a) (c | + \J a + 1 \b) {d\ 


E\ = 


-y/a |c) (a + + 1 d) (b\ 


Ei = 


a \a) (a\ + (a + 1) (|6> (6| + 



(c|) + (a + 2) \d)(d\. (7) 

This representation holds for any complex value of a. However in the following we will 
assume that a is real and positive ensuring that the Hamiltonian is Hermitian. By taking 
this representation in the expression eq. (0) we obtain an i?-matrix, denoted by R*(u), 
which satisfies the following Yang-Baxter equation 

R 12 (u - v)R* 13 (u)R* 23 (v) = R* 23 (v)R* 13 (u)R 12 (u - v) (8) 

acting on V (g> V (g> W. 
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The impurity model is constructed with generic quantum spaces represented by V and 
the impurity spaces W. To this end take some index set / = {pi,p2, ■■■■,Pi}, 1 < Pi < L 
and define 

L 

X = (g)X, 

i=l 

where 

Xi = V if i£l, 

X f - = V®W if iel. (9) 

In other words for each % G / we are coupling an impurity into the lattice which will be 
situated between the sites i and i + 1. 
Next we define the monodromy matrix 



T(u) = R 01 (u)Ro2(u)....R QL (u) 



where we have 



Roi{u) = Rai(u) for i ^ I, 

Roi(u) = RQi'(u)R* w ,(u — A) for i G I 

and we set A = a + 2. Above, the indices i' and i" refer to the two spaces in Xi (cf. eq. (|])). 
A consequence of eqs. ([3], U) is that the monodromy matrix satisfies the intertwining 
relation 

R l2 (u - v)T 13 (u)T 23 (v) = T 23 (v)T n (u)R 12 (u - v) (10) 
acting on the space V <8> V <S> A. The transfer matrix is defined by 

r(w) = tr o- T(u) (11) 

where the matrix a has entries 



a 



.i 



from which the Hamiltonian |10| is obtained through 

H = -2 -^\n(v L r(u)) 
du 



(12) 



u=0 



From eq. (pi]) it is concluded by usual arguments that the transfer matrix provides a 
set of abelian symmetries for the model and hence the Hamiltonian is integrable. In the 
next section we will solve the model by the algebraic Bethe ansatz approach. We use a 
different grading (BBF) from that adopted in [TIJ (FBB) in order to obtain the Bethe 
ansatz equations in a form which is more convenient to prove completeness of the model. 
We remark that for the case of one impurity, the equivalence of these two forms of Bethe 
ansatz equations has already been discussed in |1(J. We also present the third form (BFB) 
of the Bethe ansatz equations. 
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3. Algebraic Bethe ansatz solution 



By a suitable redefinition of the matrix elements, the operators R and R* may be written 
in terms of matrices which satisfy the Yang-Baxter equations (|3|, without Z2-grading 
(see e.g. P7j| ). These matrices read 



n a/3 



a 





























b 





c 























b 











c 











c 





b 





























a 





























b 





c 











c 











b 























c 





b 





























w 



J 

where a , /3 ( 7 and 5 ) are column ( row ) indices running from 1 to 3 and 



and 



with 



b* 

r 



1 + 2/v', 
l + 2(a + l)/u / 



2/v, 


b = 


1- 


c = - 


-2/v, 


w 


= -1 


-2/v 








1 











1 


























1 






































b* 








c* 








d 




















b* 

















e* 




















a 





























c* 








b* 








-d* 





























a 






































b* 








e* 











d* 





1 


-d* 








w* 




















-e* 

















-/* 





























-e* 








-/* 






































9* 


c 


¥ 


-2/ 




d* = 


2v^/ 


v', 


e* = 


-2Va 


+ 1/?/ 







w 



1 + 2a/v, g* = 1 + 2(a + 2)/u' 



(13) 



(14) 



(15) 



(16) 



and v' = v — a — 2. The i?-matrix ( [IB] ) acts in the tensor product of the two 3-dimensional 
spaces V ®V while i£* (|15"D acts in the tensor product of V ® W 7 , where W is the 4- 
dimensional space. Hereafter we will use these forms. 
Next we construct the monodromy matrix 



Y^u) Yfiu) Y?(u) 
Y{u) = I Y 2 \u) Y 2 2 (u) Yj(u) 
Yi(u) Yi{u) Yi(u) 



(17) 
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whose elements are operators acting in the quantum space X. The action of these elements 
is given by 

7r(Yj(u)) { $ } = (_i)(W[W]+[i][W]+W[W])3jW ( u ), (18) 
such that Y{u) satisfies the Yang-Baxter algebra without Z 2 -grading 



R 12 (u - v)Y l3 (u)Y 23 {v) = Y 23 (v)Y 13 (u)R l2 (u - v) 
The transfer matrix is expressible in terms of this representation by 



(19) 



E(-i) 



i=l 



(20) 



In the following we will omit the symbol it for ease of notation. 

For a given {a} = (a%, a 2 , a/), OLi = a,b we define the first pseudovacuum E X 

by 



,{«} 



i w 



8=1 



where 



w 



v 1 ® v aj 



for i I, 

for i = pj G I. 



An example for the case L = 10 and I = 2 is depicted below 



,{«} 



1 i ai l l l l a 2 1 1 1 1 



Setting 



a, 6 



5^>({w}) = Yf 1 Yf 2 (u 2 ) . . . . Vf* (wjv) , A = 2,3 
we look for a set of eigenstates of the transfer matrix of the form 

& = £ S^({u})v^F( M 



(21) 



where the F^ a ^ are undetermined co-efficients. We appeal to the algebraic equations 



given by eq. ( |iT| ) to determine the constraints on the variables Ui needed to force eq. (f2"T| ) 
to be an eigenstate. Although many relations occur as a result of eq. ( |i~9|) only the 
following are required: 

Y?(v)Yf(u) = a{u-v)Yf{u)Y 1 1 {v)-b{u-v)T( i {v)Y 1 1 {u) (22) 
Yf(v)Y?(u) = Y?{u)Y]\v)r%{v-u)-b{v-u)Y?{v)Y«{u) (23) 

(24) 



a{v-u)Y?(v)Y?(u) = Yf{u)Yf(v) r p a ,{v-u) 
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All of the indices in eqs. ( |2"2"| , P3] , ^) assume only the values 2 and 3. Above, the matrix 
r{u) arises as a sub matrix of R{u) and satisfies the Yang-Baxter equation for a <?/(l|l) 
invariant model. 

Using eq. (^) two types of terms arise when Y"/ is commuted through Yf. In the 
first type Y* and Y x a preserve their arguments and in the second type their arguments 
are exchanged. The first type of terms are called wanted terms because they will give a 
vector proportional to $ J , and the second type are unwanted terms (u.t.). We find that 



N 



Y?(v)& = a(v) L J] a{ Ul - v)& + u.t. 
i=i 

where we have used the fact that 

Yl(v)v {a} = a(v) L v {a} 
For the case of the previous example we have 

11 OL\ 1111 « 2 1111 

Y l v{a} 



(25) 



(26) 



a, b 



11 a x lili a 2 1111 
Similarly, for i = 2, 3 we have from eq. (|23"D (no sum on z) 

= ^'>(M)^(^4 { |>(^M)^>^ jQ} + u.t. 
= ^(M)^^, M)^' } i^ } + u.t. 



where 
Above 



fc{/3} 



*(«,{«}) = r i(u-Ui)ro2(i;-U2)...rojv(u-wjv), 



r OlW r 02W-- r oA^J 



and r*(w) is the g/(l|l) i?-matrix defined from R*(u) by restricting the indices to a, 7 = 
2, 3 and /3, 5 = a, b in (^) in analogy to r(u). It is important to observe that the space 
X is closed under the action of the elements YJ(u), i,j = 2,3 which generate a <jrZ(l|l)- 
Yang-Baxter algebra. In fact 

Yj(u)vM = t$$ } (u)vW 

which follows from the fact that the YJ(u) i,j = 2, 3 act trivially on the vector v 1 in the 
sense 



Y^iuW = Y?(u)v 



Y^v 1 = Y 2 3 {uy = 
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An example of the action of Y 2 (u) and Y£{u) on the first pseudovacuum is illustrated 
below 



11 a x 1111 a 2 1111 



Y 2 v {a} 



11 a x 1111 a 2 1111 



a.i = a,b 



11 ax llll a 2 llll 



y 3 M a > 



11 ai llll a 2 llll 
The contribution to the eigenvalues of the transfer matrix is 

3 



«j = a, 6 



i=2 



(27) 

At this point we need to perform a second-level, or nested Bethe ansatz procedure to 
diagonalize the matrix 



i]+[i][{/3,a}]Tj{/3,a} 



i=2 



where we have used the fact that F{p a \ = unless [j] = [{(3, a}]. The above matrix 
is simply the transfer matrix for a g/(l|l) invariant system acting in the tensor product 
representation of N copies of the vector representation with inhomogeneities {u} and I 
copies of the impurity sub-representation. 

We diagonalize Tx(v) in analogy to t(v) ( p0[ ) by constructing the nested monodromy 
matrix 

yl{u) yl(u) 

.yi(u) yl{u) 

satisfying 

r 12 (it - v)yx3(u)y 23 (v) = y 23 (v)yx3(u)r 12 (u - v). 



y{u) 



:28) 



From the above set of relations we will need the following 

yl( v )yl( u ) = a ( u - v )yl{ u )yl{ v ) - K u - v)y 2 (v)y 2 (u), 
yl{ v )yl( u ) = ~ a ( u - v )y2( u )yl( v ) - K v - u )yl{ v )yl{ u )i 



y 3 2 (v)y 3 2 (u) 



—a(u — v) 



yl(u)y 3 2 (v). 



a(v — u) 

Proceeding similarly as before, we look for eigenstates of the form 

<t> = yUix)y 2 (i2)-ylhM)0 



(29) 



(30) 
(31) 

(32) 
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with the second-level pseudovacuum 9 given by 

6 = S {2 \{u})v {a \ 

It is an eigenstate of y\ and y\ 



N 



y%(v)9 = n <v - Ui)d , yl(v)6 = [w*(v)] l 9 . 



(33) 



i=l 



Using (|30|j3l|) it follows that 



T\(y)(j) = Ai(v)<f> + u.t. 



with 



N M M 

A l( v ) = II a ( V - U i) II - V) - [W*(v)f Y[ a{lk - V) 



i=l k=l k=l 

The unwanted terms cancel provided the parameters satisfy the Bethe ansatz equations 
(BAE) 

N 

J] a(j k - Ui ) = [w*( lk )f, k = 1, 2, .., M. 



(34) 



Combining this result with eq. (^) we obtain for the eigenvalues of the transfer matrix 
eq. (P) 



N 



A(v) = a(v) L J[ a(ui - v) + Ai(u). 



(35) 



i=l 



Cancellation of the unwanted terms in (|25 , 27D leads to a second set of BAE which are 



t „ a{ Ui - u h ) 
i{u h ) 11 



M 



tA a(u h - Ui ) k=1 
Using ( |I^JT6| ) and making the change of variables 

u h — > iu h + 1 
7fc -> ilk + 2 



IIa(7fc-Wfc) & = 1,2,...., A. 



(36) 



we obtain 



f u h + i 
\u h - % 



» [ u h -u k + 2i^u h ->y s -i h = 1 N = L + 1 _ 
k -t h u h -u k - 2i f* x Uh-js + i 



-q Is ~ u h + i (is- ioc\ l 



S = l,...M = L + 2-ni i -n l 



(37) 



(3* 



£=iis-Uh-i \is + i<yj 

In the absence of impurities (limit I — > 0) we recover the form of the BAE first derived 
by Sutherland || and later by Sarkar for the usual i — J model. The basic procedure to 
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solve ( P?| , PBD is to adopt the string-conjecture, which means that the Uh appear as strings 
and all roots 75 are real 

<j = ul + i(n + l-2j) ;j = l,...,n a = l,...,N n n = 1,2,3,... 
Js = real 5= 1,2,3, ...M (39) 

where is the position of the center of the string on the real w-axis. The number of 
n-strings N n satisfy the relation 



Applying the string conjecture (^) in fl3~7|) and (|38D and taking the product over j in (|37|) 
and (133) we obtain 



e 



L '"^ -UUEnU<-uf)f[e^(^-^) (40) 

m j3 S 

nn4^)=e-(?) (4i) 



n 



where e(x) = —. and 



X \ 


» 2 


( X 


n — m 






n — m 



E nm {x) = e ■ r e 1 — - ) ...e — ■ e — ■ ) , n^m 

\n + m — 2 / \n + m; 



^fMlJ-^feMsJ- "= m (42) 

Taking the logarithm of these equation and using ln(e (x)) = ? arctan (x) — In (— 1) we get 

^ Ur) - E E Q-m (< - «?) + E ' ( = 2 - J « (43) 



" ' m /3=1 5=1 ^ n 



53 £ e (7*_J^\ + w( 75 } = 2?r Js (44) 



where 8(x) = 2arctan(x) and 



Qnm(x) = 6 U^--)+29 X -—- + +2e( . * Wflf— g— ), „ 



|n — m|/ \ |n — m\ + 2) \n + m — 2J \n + m 

■2H (?) +2 e(-)+... + 2e( +e(^-), n = m 



2 J V4/ \2n-2j \2n 



Hence the solutions of (|37|) and (|3~8| ) are parametrized in term of the numbers 7™ and 
J/3. The numbers J™ and J/3 are limited to the intervals 



\a < c = ^ ( L + M - E wvm - 1) 



(46) 
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\Js\ < J max — - N n + l — 2^ 



(47) 



where t nm = 2min(n, m) — 5 n m . 

As was shown in the papers ||, ||| two other forms of the Bethe ansatz exist which are 
obtained by choosing a different grading for the indices of the gl(2\l) generators. Recall 
that the above calculations were performed with the choice 

[1] = [2] = 0, [3] = 1. 

Choosing 

[1] = 1, [2] = [3] = 

yields the eigenvalue expression 

N 



a(v) = -HKtO]V(tO] , IM«-«*) 

, N M M 

+ [b*(v)f (l[a(v-u i )l[a( lk -v) + l[a(v- lk )) (48 



i=l k=l k=l 



subject to the BAE 



[- w ( u i)] ( 777 \ ) = Y[a(lk-u t ), i = l,2,...,N 

V \ u i) / fc =1 



i H 7fe 7fc | = /i=i,2,...,m. 



In the limit / — > we recover Lai's form of the BAE 0] (see also 0). Alternatively, 
choosing 

[1] = [3] = 0, [2] = 1 

yields the eigenvalue expression 



N M 



A(v) = a(v) L Y[a(u l -v) + [b*(v)] l l[a(v- lk ) 

i=l k=l 
N M 

- lf*(v)} l I[a(u i -v)l[a(v- lk ) 



with the BAE 



a(Ui 



i=l k=l 



M 

l L = [/»]' II a (^-7fc) i' = l,2,...,JV, 



(^)'ft«(«*-7 fc ) fc = l,2,...,M. (49) 



fc=i 

l N 
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Finally, from the definition of the Hamiltonian eq. (12) we see that the energies are 
given by 



dv 

Using the eigenvalue expression eq. (B3) we obtain 



v=0 



N 



E 



+ < 



where the Ui are solutions to the equations ( p7| , p8D . 



4. Highest weight property 

Next we wish to show that the eigenstates constructed in the previous section are in fact 
highest weight states with respect to the underlying supersymmetry algebra gl{2\l). The 
highest weight property of the Bethe states has been proved for many models, such as the 
Heisenberg chain [23] and its generalized version pB|, the Kondo model [p6 |, the usual 
t — J model 0, the Hubbard chain [|4| and its gl(2/2) extension [25|. In all of these 
examples the Bethe ansatz was performed using a unique reference state. However it was 
shown recently that the highest weight property can also hold for cases where the Bethe 
ansatz solution is obtained from a subspace of reference states [|TB . 
Let us begin by considering 



£ E 2 3 S^({u})v^F : 
{/»,«} 



•3 



By means of the nesting procedure we know that the co-efficients F^ , are such that we 
have the following identification of states 



for a suitable solution of the BAE. By comparing eqs. (|6|, |T8| , |29D it is possible to determine 
algebraic relations between the elements of the Yangian algebra and the supersymmetry 
algebra. For our purposes we need the following 



[El yl{u)V p = -y 2 a (u)2 + y!(uT ( 
Noting that E%w = it is evident that we may write 



(50) 



M 



E lyl{ii)--yl{ , iM)w = J2 x h x h 



h=l 



with 



x h = 2/2(71) yl{ih~i)yl{ih+i)--yl{iM)w 

and the Xh some yet to be determined co-efficients. To find Xh we write 



yl{ii)--vl{iM)w = n 



yl{ih)x h 
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where we have used eq. Now by using the relations ([J(],[n],^t]) and looking only for 

those terms which give a vector proportional to Xh we find that 



J=i a(7? -7/0 



M N M 

i w *bh)] 1 n - 7^) - n a (7ft - u d n °(7* - to 

k^h i=l fc^/i 



which vanishes because of eq. (|34|) . Thus we see that 

£f $ j = 0. 

Next we consider the action of E\ on . Using eqs. (|6| JT9|JT8|) we find the commutation 
relation 



Y{*(u)] — 6%Y{(u) — Yr 



[u 



(51) 



As before, since we can write the general expression 

E\& =J2z hjP Z h>f3 



where 



Z hiP = SK>({ui)})SWH{ut)})vl"W VM 
and for any vector {w} we have 



{ Wh } = ( Wlj yj 2l lW A _i), {w£} = {W h+1 , 

To calculate z^p we begin by writing 

& = ^ } ({%})^K)^ +} (K})^ W ^ 



W7V 



/i.-l 



I] a(«i - u^-^^C-Tifc, {-«;: })Fx 7 K)SW(K })^ +} ({<})^ W ^} 



i=i 



where we have used the relation eq. (24). Now applying eq. (|5l|) and using the relations 
(fl^,^) to determine the terms which give a vector proportional to Zh,p we find that 

(N N M 

a{uh) L n a ( u i ~ - n ~ n a (7fc - 
j^ft i^/i fe=l 



which vanishes as a result of eq. (0). We then conclude that 

E\& = 0. 

Finally, using the fact that 



it is deduced that 



E l = [ E 2i E V\ 



El& = 
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which completes the proof that the Bethe states are gl(2\l) highest weight states. We 
observe that this property can also be proved for the other two choices of gradings in a 
similar way. 

By acting with the gl(2\l) lowering operators on the Bethe states additional states 
are obtained which are necessarily eigenstates of the Hamiltonian as a consequence of 
the gl(2\l) invariance. Before proceeding to counting the number of eigenstaes, it is 
important to assert that each of the Bethe states belongs to an irreducible submodule 
of the full tensor product space. Recall that unlike the case of Lie algebras, there exist 
representations of Lie superalgebras which can be indecomposable but not irreducible |29| . 
However this situation does not occur in the decomposition of the tensor representation 
here, which follows from the fact that both the fundamental representation and the family 
of four dimensional representations with a > belong to the class of type I unitary (or 



star) representations defined in []29| and classified in P0 |. A general theorem states P?J 
that tensor products of type I unitary representations are always completely reducible 
into representations which are also of type I. Moreover, since the module which acts for 



the impurity spaces is typical it follows from [^T| (see Proposition 2) that all the modules 



occuring in the decomposition of the tensor product space are in fact typical. 



5. Completeness of the Bethe states 

In this section we show how to construct a complete set of eigenvectors for the t — J 
hamiltonian for an arbitary chain of length L with Z-impurities. This is obtained by 
combining the Bethe ansatz equations with the supersymmetry of the model. 

From the previous section we know that all collections J$} where the J's and J's 
are pairwise different specify all the Bethe vectors. The number of admissible values for 
the I"'s and the J^'s (for fixed {N n } and M) is (2J" ax + 1) and (2J max + 1), respectively. 
Taking into account that many different string configurations N n give the same number 
of roots N , the number of possible Bethe vectors for fixed N, M is given by 



z(n,m) = f 2Jm r; +1 )n 

{N n } \ m In 



max 1 



(52) 



where the sum over {N n } is constrained to ^Z n nN n = N. It is convenient to introduce 
the quantity q = J2 n N n - Then we can write this sum as follows 

m, m) = £ (" + < - x ) £ n ( L ~ E ™ Y" + M ) • (53) 

g=0 V iV1 J {N n } n \ iV « / 

where the inner sum is constrained to fixed values of N and q. This expression resembles 



the one calculated by Bethe in the isotropic Heisenberg model |22|, ^3|] and can be simplified 
to 

Z(N, M) = + M — 2N + I fq + I — l\ fL + M — N + l\ (N — l\ 
1 J ^ L + M-N + l\ M )\ q J\q-l) 1 ^ 
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The total number of Bethe vectors is obtained by summing Z(N, M) over all N, M. 
From the fact that the Bethe vectors are highest weight vectors, by acting with the 
gl(2\l) lowering operators EJ- (i > j) on the Bethe states we obtain additional states. Each 
Bethe state (with fixed N, M) is the highest weight vector in an irreducible multiplet of 
dimension [ f29fl 

d(N, M) = 8(S Z + 1/2) = 4(L — 2N + M + 1) 

Notice that here, in constrast to the usual t — J model || only typical gl(2\l) multiplets 
occur as discussed in the previous section. With these considerations, the total number 
of eigenvectors is 



L+M 

L+21 ~+~ 



Z =H E d(N,M)Z(N,M). 

M=0N=M-l 

Above the ranges of M and N are obtained by demanding that the magnetization S z = 
| (n-j- — ii i) = | (L — 2N + M) and the number of electron Q = n-f + nj_ = L + 21 — M are 
restricted to intervals < S z < ^ and < Q < L + 21. 



L + M 



f ± HL-2N + M+l) L + M - 2M + 1 x 
kUnk-i ' L + M-N + l 

U 9+ L-Y +M ; N+1 )(::t) 



q=0 



Using the identity 



. M ) f^A k \M-k 



k=0 



eq.(E3) can be written as 



L+M 

(I - i\ L +^ 1 a L + M - 2M + 1 

z = y [ , ) y y a (l - 2N + m + 1) — — x 

»( q \(L + M-N + l\(N-l\ 

§ M-J a L-l • (56) 



g=0 



Employing some combinatorics, we find 

L+M 

l-l // _ l\ £+2* — T - 

^ = E . E E 4(L-2iV + M + l)x 

fc=0 \ K ) M=0N=M-l 

L + M -2N + \(L + M -N + l\f L + k \ 
L + M-N + l \ M-k )\N -M + kj 

After rearrangement this expression can be rewritten as 

L + M 

l-l (] _ \\ L+21 

Z = E . E E HL-2N + M+l)x 

k=0 \ / M=0N=M-l 

UL + M -N + l\( L + k \ ( N \(L + k^ 
IV M-k J \N-M + k) ~ \M-kj \N - 1 



(57) 



(5f 
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Making the substitution iV — > x + M — I we obtain 



L+2Z — 2— 



ST/! E E 4(L-M + 2/-2o; + 1)x 



fc=0 



M=0 z=0 



\(L + l-x + l\( L + k \ (x + M-l\( L + k 
{{ M-k )\x -l + k) ~ \ M-k )\x + M -l-l y 

which after some manipulation turns out to be 

l-l /] _ 1 \ L+21 

z = 4 E( k J E (m-* + i)x 



, (59) 



fc=0 

f L — M 



M=0 



+1 



-+l 



- " / L + l-x + l\f L + k \ (x + M L + k 



E 

z=0 
2-1 

-4E 



M — A; + 1 / Vx - I + k 



+ E 

/ rr — n 



x=0 



Jtf-ife + lHx + M-Z-l, 



7- r 



L+22 



E (£ + *)x 



k=0 

L-M 



M=0 



+1 



E 

x=0 



L + l — x + i\f L + k -I 



-+i 



M-k 



,x - I + k - 1, 



E 

z=0 



' L + k — I \(x + M-f 
x + M -I- 1/ I M-k , 



(60) 



Substituting x ^ L — x — M + 2/ + 1 in the second and fourth terms of the equation above 
we get 



i-l // l\ L+2Z L-M+22+1 /r ~ i 7 i 1 



k 

k=0 \ h J M=0 
l-l _ ]\ L+21 



x=0 



-4E 

fc=0 



^ . ,, L - M + 2l+1 [L-x + l + 1 

E ( L + k ) E 



M=0 



x=0 



M-k + 1 ){x-l + k) 

L + k — 1 N 
M-k ){x-l + k-l t 



(61) 



We now find 



froV k JjfcoKM-kJdP 

t^a-iy+? (L + k\_d_ 

and using the binomial formula we obtain 



L-M+2k 

E 

x=0 



'L-M + 2fc\ 
x — I + k J 



J p=l 



L -g- 2 YL-M + 2^ p:c _ L _ z _ 1 



x=0 



-l+k-1 



(62) 



p=i 



2=4 

dP 



^ 1 Vi-p- 1 )p-*- i E J(i + ^) L - M+2fc 



E 



A: 



Lfe=0 \ / M=0 

Above the sum over M can be performed giving 



M-k, 



(63) 



j p=i 



Z = 4— 
dP 



i-i 



(i-^)E 



fc=0 



7-l N 



p- k ~ L (2 + P) 



L+k 



(64) 



p=i 
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Now the sum over k can be evaluated yielding 



4 A 

dP 

which finally leads to 



P-')P 



-L-l+l 



(2 + P) L {2 + 2P) 1 



p=i 



3 L A l . 



(65) 



(66) 



Thus we have shown that the number of eigenvectors of the t — J hamiltonian with / 
impurities is 3 L 4 l , which is precisely the number of states in the Hilbert space of a chain 
of length L with / impurities. We point out that this is the first time that completeness 
of the Bethe vectors has been demonstrated for an impurity electronic model. 



Acknowledgements 



JL thanks the Fundagao de Amparo a Pesquisa do Estado do Rio Grande do Sul and 
Australian Research Council for financial support. He also thanks the Instituto de Fisica 
da UFRGS for their kind hospitality. AF and APT thank CNPq-Conselho Nacional de 
Desenvolvimento Cientffico e Tecnologico for financial support. 

References 

[1] N. Andrei and H. Johannesson, Phys. Lett. A100 (1984) 108 

[2] P. Schmitteckert, P. Schwab and U. Eckern, Europhys. Lett. 30 (1995) 543 

[3] H-P. Eckle, A. Punnoose and R. Romer, Europhys. Lett. 39 (1997) 293 

[4] C.K.Lai, J. Math. Phys 15 (1974) 1675 

[5] B. Sutherland, Phys. Rev. B12 (1975) 3795 

[6] P. Schlottmann, Phys. Rev. B36 (1987) 5177 

[7] S.Sarkar, J. Phys. A24 (1991) 1137 

[8] F.H.L. Essler and V.E. Korepin, Phys. Rev. B46 (1992) 9147 

[9] A. Foerster and M. Karowski, Phys. Rev. B46 (1992) 9234; Nucl. Phys. B396 (1993) 
611 

[10] G. Bediirftig, F.H.L. Essler and H. Frahm, Phys. Rev. Lett. 77 (1996) 5098; Nucl. 
Phys. B498 (1997) 697 

[11] A.J. Bracken, M.D. Gould, J. Links and Y.-Z. Zhang, Phys. Rev. Lett. 74 (1995) 
2768 

[12] J. Abad and M. Ri'os, Phys. Rev. B53 (1996) 14000 



19 



[13] J. Abad and M. Rios, J. Phys. A30 (1997) 5887 

[14] M.P. Pfannmuller and H. Frahm, J. Phys. A30 (1997) L543 

[15] J. Abad and M. Rios, "Exact solution of a Electron System Combining Two Different 
t - J Models" , |cond-mat/9806106| 



[16] J. Links and A. Foerster, J. Phys. A32 (1999) 147 

[17] M. J. Martins, "Integrable mixed vertex models from braid-monoid algebra" , preprint 
UFSCAR-TH-98-12 



[18] P.-A. Bares, Exact results for a one dimensional t — J model with impurities |cond- 
| mat/ 9412011 



[19] P. Schlottmann and A. A. Zvyagin, Phys. Rev. B55 (1997) 5027; B56 (1997) 13989 

[20] Y. Wang, J.-H. Dai, Z.-N. Hu and F.-C. Pu, Phys. Rev. Lett. 79 (1997) 1901; 
Z.-N. Hu, F.-C. Pu and Y. Wang, J. Phys. A31 (1998) 5241. 

[21] H.-Q. Zhou and M. D. Gould, Algebraic Bethe ansatz for integrable Kondo impuri- 



ties in the one-dimensional supersymmetric t — J model, cond-mat/9809055 
H.-Q. Zhou, X.-Yu Ge, J. Links and M. D. Gould, Graded reflection equations al- 
gebras and integrable Kondo impurities in the one-dimensional t — J model, |cond- 
mat/9809056; . 



[22] H. Bethe, Z. Phys. 71 (1931) 205 

[23] L. D. Faddeev and L. Takhtajan, Zap. Nauch. Semin. LOMI 109 (1981) 134 

[24] F.Essler, V.Korepin and K. Schoutens, Phys. Rev. Lett. 67 (1991) 3848; Nucl. Phys. 
B372 (1992) 559; Nucl. Phys. B384 (1992) 431 

[25] K. Schoutens, Nucl. Phys. B413 (1994) 675 

[26] A.N. Kirillov, J. Sov. Math. 34 (1985) 2298 

[27] G.W.Delius, M.D.Gould, J.R.Links and Y.-Z. Zhang, Int. J. Mod. Phys. A10 (1995) 
3259 

[28] A.N. Kirillov, J. Sov. Math. 36 (1987) 115 

[29] M. Scheunert, W. Nahm and V. Rittenberg, J. Math. Phys. 18 (1977) 146. 
[30] M.D. Gould and R.B. Zhang, J. Math. Phys. 31 (1990) 1524 
[31] M.D. Gould and J. Links, J. Phys. A30 (1997) 1613. 



20 



